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Abstract. We examine square and kagome artificial spin ice for colloids confined in 
arrays of double-well traps. Unlike magnetic artificial spin ices, colloidal and vortex 
artificial spin ice realizations allow creation of doping sites through double occupation 
of individual traps. We find that doping square and kagome ice geometries produces 
opposite effects. For square ice, doping creates local excitations in the ground state 
configuration that produce a local melting effect as the temperature is raised. In 
contrast, the kagome ice ground state can absorb the doping charge without generating 
non-ground-state excitations, while at elevated temperatures the hopping of individual 
colloids is suppressed near the doping sites. These results indicate that in the square 
ice, doping adds degeneracy to the ordered ground state and creates local weak spots, 
while in the kagome ice, which has a highly degenerate ground state, doping locally 
decreases the degeneracy and creates local hard regions. 
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Artificial spin ice systems constructed with nanomagnetic arrays m El a n El El m E], 
vortices in nanostructured superconductors 0HQ1III], or soft matter systems mm 
fl4l M El H?] have attracted growing attention as outstanding model systems in which 
different types of ordered and degenerate ground states can be realized ms, as well as 
various types of avalanche dynamics [5j d8l [19], return point memory [20], and a variety 
of thermal effects [2_1, 22ll23ll2Tl 251126] , Another key feature is that many artihcial spin 
ice systems are constructed on size scales at which the microscopic degrees of freedom 
can be imaged directly [8j. These systems are called artihcial ices since their ground 
state can obey what are called the ice rules that were first studied in the context of a 
particular phase of water ice. Here, on each bond between oxygen vertices, two protons 
are localized close to the vertex and two are localized far from the vertex, creating what 
is called the “two-out, two-in” rule [27]. Since there are many ways to arrange the 
effective charges or spins while satisfying the ice rule constraint, the ground state is 
highly degenerate so that even at T = 0 there is an extensive entropy. This picture has 
also been applied to certain classes of atomic spin materials with pyrochlore structures, 
known as spin ice systems [23 E3 SO, El]. 

In the artihcial spin ices, a set of geometric constraints is imposed on the system 
via the arrangement of the nanomagnetic islands [8] or arrays of double-well traps for 
vortices [3 Es .H] or optical traps for colloids [121 El • Here, each trap or nanoisland 
plays the role of an effective macroscopic spin. For nanomagnetic arrays, the spin 
direction is defined to point in the direction of the magnetic moment, while for the 
colloidal or vortex systems with double-well traps, the spin is defined to point toward 
the end of the well that contains a particle. The islands or traps are arranged in a 
square geometry as shown in Fig. [T] to create artihcial square ice, or in a hexagonal 
arrangement to realize artihcial kagome ice. In the square ice arrangement there is an 
ordered ground state in which each vertex obeys the two-in, two-out ice rule. The lowest 
energy excitations in this system take the form of vertices with three spins in or three 
spins out to form a monopole with charge —1 or +1, while higher energy excitations are 
vertices with four spins in or four spins out, forming monopoles with charge +2 or —2 
HI IS]. There are also ice rule obeying states known as biased states that have somewhat 
higher energy than the ground state, and often pairs of monopoles can be connected 
by a string of biased state vertices. The kagome ice ground state is not ordered but 
does obey the ice rules, which in this case are two-in, one-out or one-in, two out. Here, 
monopoles consist of vertices with three spins in or three spins out. There are additional 
artihcial spin ice geometries in which the introduction of different constraints or rules 
produces various types of ordered or quasiordered ground states [32[ [33, Mj • 

It is also interesting to consider how quenched disorder can affect the ice states. 
In the magnetic spin ice systems, disorder can arise as a dispersion in the energy of 
the barrier that must be overcome to hip an individual effective spin. It should also 
be possible to create positional disorder in the system by shifting the nanoislands away 
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from their regular lattice positions. In square ice systems, disorder can lead to the 
formation of domain walls, with non-ground state vertices dividing the system into 
separate regions of different ground state vertices PEI 1351136]. Additionally, defects on 
individual nanoislands can affect the interactions between monopoles [37] . One method 
for characterizing the effects of disorder is by applying a cyclic held or drive sweep to 
generate hysteresis loops between two biased ground states. Libal et al. constructed 
hysteresis loops for colloidal artificial spin ice containing disorder in the heights of the 
barriers the colloids must overcome to hop from one side of the trap to the other [ffl- 
They observed that the square ice forms domain walls that gradually coarsen during 
repeated cycling of the drive, increasing the fraction of the sample that is in a ground 
state configuration. In contrast, the kagorne ice contained no domain walls, and the 
number of non-ice rule obeying vertices at zero bias was almost constant, indicating 
that the system underwent little to no coarsening. In the square ice the dominant mode 
of defect motion was propagating domain walls, while in the kagorne ice it was individual 
hopping of defects, which became pinned and stationary rather than annihilating ED). 

For superconducting or colloidal artificial ices, it is possible to add disorder in 
the form of an effective doping by placing doubly occupied or unoccupied traps in the 
sample. A doubly occupied trap contains an effective spin that points toward both 
vertices on either end of the trap at the same time, while a doubly unoccupied (empty) 
trap contains an effective spin that points away from both vertices on either end of the 
trap at the same time. In magnetic spin ice systems such configurations are not possible; 
however, for water ice, additional protons can be added or removed at individual bonds, 
allowing for similar arrangements. In atomic spin ice systems it is possible to create 
what are called stuffed spin ices by chemical alterations, against which the spin ice rules 
have been shown to be robust [38] • There are also studies on diluting spin ice systems 
where the dilution can generate new degrees of freedom [39] . 

In this work we consider colloidal spin ice for both square and kagorne geometries 
that have been doped by adding extra colloids to create doubly occupied traps. We 
examine how the ices disorder under the application of thermal fluctuations. Colloidal 
systems are ideal for studying such doping effects since large scale arrays of optical traps 
can be realized and the number of colloids per trap can be precisely tuned [40, 41]. An 
example of the square ice geometry appears in Fig. 1, which illustrates the local energy 
configuration generated by the double-well traps, each of which captures a single colloid 
except for the doped site at the center of the image which contains two colloids. The 
doped site creates a geometrically necessary three-in, one-out higher energy vertex. 

2. Simulation Method 

We perform two-dimensional Brownian dynamics simulations of colloids in double-well 
traps using the same techniques applied in previous work on colloidal spin ice systems. 
Our square ice sample contains N = 4232 double-well traps arranged in a 46 x 46 lattice 
to form 46 x 46 = 2116 vertices, while our kagorne ice contains N = 4200 double-well 
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Figure 1 . Illustration of a doping defect in colloidal artificial square ice. Each trap 
confines single colloid except for the center trap, which is doped and contains one 
colloid at each end. The energy landscape created by this configuration is plotted in 
the lower plane of the figure. 


traps arranged in a 20 x 35 lattice to form 20 x 35 x 3 = 2100 vertices. The simulation cell 
is of size 90ao x 90ao for the square ice and 120a 0 x 121.24ao for the kagome ice system, 
where ao is the simulation unit of distance, and it has periodic boundary conditions in 
both the x and the y directions. The elongated traps are 1.8a 0 long and 0.4a 0 wide 
and are placed a center-to-center distance a = 2.0ao from each other so that they never 
overlap. The overdamped equation of motion for colloid i is: 


V 


dRj 

dt 


pcc , yl 
i ' i 


+ Ff 


( 1 ) 


where the damping constant rj = 1.0. The colloid-colloid interaction force has 
a Yukawa or screened Coulomb form, F^ c = — F 0 q 2 Y^j ViY(r^) with Vfaj) = 
(l/ r u) exp(—rjj/r 0 )fjj. Here r tj = |iy - r,-|, % = (r* - rr i(i) is the position 
of colloid i(j ), F 0 = Z* 2 /(47ree 0 ), Z* is the unit of charge, e is the solvent dielectric 
constant, q is the dimensionless colloid charge, and r 0 is the screening length, where 
r 0 = 4a 0 so that interactions extend as far as second neighbors of each trap. We 
neglect hydrodynamic interactions between colloids, which is a reasonable assumption 
for charged colloids confined in traps that remain in the low volume fraction limit at all 
times. The thermal force F T is modeled as random Langevin kicks with the properties 
(Fj) = 0 and (Ff (t)Fj(t’)) = 2 rjk B TSijS(t—t'). We heat the system by slowly increasing 
T in increments of 5T = 0.05 steps from T = 0 to T = 18. At each temperature we 
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allow the system to equilibrate for 50, 000 simulation time steps. As the temperature 
increases, the colloids begin to hop between the two minima of the double-well traps. 

3. Square Ice: Local Screening and Melting 

In Fig. [2](a-c) we highlight a portion of the square ice system with a doping ratio of 
x = 0.0236, where 2.36% of the traps are doubly occupied. The doped colloids are red 
and the different vertex types are indicated by the colors shown at the top of the figure. 
In the absence of doping, the fraction Pqs of vertices in the ground state at T = 0.0 is 1.0, 
and the system becomes disordered with Pgs — 0.75 near T m 14, as shown in Fig. |3](a). 
We use the nomenclature N2 gs for the ground state vertices; N2bi ase d for the biased 
two-in, two-out states; N3 for three-in, one-out -1 monopole vertices; N1 for three-out, 
one-in +1 monopole vertices; N4 for a four-in -2 double monopole; and NO for a four- 
out +2 double monopole. Figure 2(a) illustrates a low temperature state with T — 1, 
where each doubly occupied site is screened by an adjacent N3 vertex but there are no 
N1 excitations in the sample. This is in contrast with the zero doping limit, where for 
each N3 monopole excitation there must exist a compensating N1 monopole excitation. 
Such ±1 excitation pairs can be created through a single spin flip, while subsequent 
spin flips make it possible for the magnetic charges to move some distance away from 
each other through the lattice. An emergent attractive force will arise between the two 
opposite monopoles that can be described by a modified Coulomb law for magnetic 
charges, and a Dirac string of N2 biased vertices that connects the monopoles becomes 
more energetically costly as the monopoles move apart and the string becomes longer. In 
the doped system, introducing a double defect creates an isolated N3 monopole that does 
not have a compensating N1 monopole. This indicates that the square spin ice responds 
to the the doping by locally screening the extra added charge through formation of — 1 
excitations next to the double defects, as shown in Fig E^a) , where every double defect 
has a screening N3 defect next to it. 

In Fig. [2](b) we show the same x = 0.0236 square ice system at T = 5.0. The 
background ordered ground state remains frozen, but two types of additional excitations 
emerge near the doubly occupied sites. In the first excitation, the screening N3 monopole 
starts to move away from the doped site but remains connected to it by a string of 
N2bi a sed vertices, as shown in the upper left portion of Fig. [2](b) where an N3 defect is 
separated from the doped trap by a single N2 b i ase d vertex. In the second excitation, a 
string of ±1 monopoles is created, as shown in the center right portion of Fig. [2] where 
the single N3 defect has turned into a pair of N3 defects separated by an N1 defect. In 
Fig. [2](c) , which shows the same doped square ice sample at T = 12, the ground state 
regions away from the doubly occupied traps are still ordered; however, a larger number 
of strings of N2 b i aS ed vertices or N1-N3 defect strings have formed near the doped sites. 
This result shows that the doping introduces extra topological charges that can serve 
as nucleation sites for the thermal wandering of monopoles at temperatures well below 
the bulk melting temperature. The doping sites can be viewed as local weak spots that 
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Figure 2. Images of small portions of the square and kagome ice systems showing the 
double-well traps (open ellipses), colloids in singly-occupied traps (filled black circles), 
and colloids in doubly-occupied traps (filled red circles). The larger colored circles 
indicate the vertex type as illustrated at the top of the figure: for the square ice, NO 
(blue), N1 (green), N2bi aS ed (yellow), ground state N2 gs (white), N3 (orange), and N4 
(red); for the kagome ice, NO (blue), ground states Nl gs and N2 gs (white), and N3 
(orange). (a,b,c) The square ice system at a doping of a = 0.0236. (a) At T = 1.0, 
below melting, each doped site is screened by an N3 monopole, (b) At T = 5.0, some 
thermal wandering of the N3 sites can occur, and N1 states can form at the doped 
sites, (c) At T = 12 the regions away from the doped sites remain in the ground state 
while local melting occurs at and near the doped sites. (e,f,g) The kagome ice system 
at a doping of x = 0.095. (e) At T = 1.0, below melting, the ground state absorbs the 
doping charge without forming defects by increasing the ratio of N2 gs to Nl gs vertices, 
(f) At T = 12.0, N3 monopoles form in regions away from the doped sites, (g) At 
T = 15.0, N1 monopoles begin to appear. 
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Figure 3. (a) The fraction of ground state vertices Pgs vs T for square ice samples 

at doping fractions of x = 0, 0.00236, 0.00472, 0.00945, 0.0709, 0.0945, 0.1182, 0.142, 
0.189, 0.2363, 0.2835, 0.33, and 0.472, from top to bottom. A two step disordering 
process occurs, (b) Pqs vs doping fraction x at T = 0 (blue circles) and T = 2.0 (red 
squares). Blue dashed line: a fit to Pgs = 1 —a:. Red dashed line: a fit to Pgs = 1 — 2x. 


induce local melting. We also observe similar behavior at higher doping levels in the 
regions that are not adjacent to the doping sites. 

In Fig. [3(a) we plot Pgs, the fraction of vertices that are in the ground state, versus 
T for the square ice system from Fig. [2] for doping fractions x ranging from 0.0 to 0.472. 
For a given temperature, Pqs decreases approximately linearly with increasing x as each 
doped site is screened by the formation of a non-ground state vertex. We observe a two- 
step disordering process for finite doping, as indicated by the drop in Pqs near T = 1.5, 
which becomes more pronounced as x increases. At around this temperature of T — 1.5, 
the N3 screening monopoles begin to hop between neighboring sites as the individual spin 
degrees of freedom within the vertex begin to flip, and as they hop, they create additional 
non-ground state vertices in the form of N2bi aS ed or N1 monopole sites, depressing the 
value of Pgs- As T increases further, additional thermally-activated defects appear as 
the system approaches the clean melting temperature, and the resulting decrease in Pqs 
begins at lower values of T as the doping level x increases due to interactions between 
vertices surrounding the randomly placed doping sites. In Fig. [3(b) we plot Pqs versus 
x for T = 0 and T — 2. At T — 0, if each doped site created a single non-ground state 
vertex, the curve would follow the blue dashed line, which is a fit to Pqs = 1 — x. Instead 
the curve drops below this value, indicating that as the doping fraction increases and 
some vertices begin to interact with more than one doped trap, more than one screening 
vertex can form per doping site. At T — 2, just above the temperature at which the 
screening N3 monopoles become able to hop to an adjacent site, there is on average one 
additional defected vertex generated for every screening vertex in order to permit this 
hopping, as indicated by the red dashed line, which is a fit to Pgs = 1 — 2x. 
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Figure 4. Pqs vs T for the kagome ice system from Fig. [2] at x = 0.0, 0.095, 0.0238, 
0.0476, 0.0714, 0.0952, 0.119, 0.143, 0.19, 0.238, 0.286, 0.33, and 0.476, from top to 
bottom, (b) Pqs vs x at T = 0.0 for the kagome ice (green circles) and square ice (blue 
squares), showing that doping only weakly affects the ground state configurations of 
the kagome ice. (c) Pqs vs x for the kagome ice at T = 17, with error bars of width 
0.001. At this temperature, low levels of doping can reduce the defect hopping and 
diminish the disorder in the ground state configuration compared to the undoped case, 
(d) Pqs vs x for the square ice at T = 17, with error bars that are smaller than the 
symbols. Addition of doping monotonically decreases the ground state order in the 
system. 


4. Doped Kagome Ice States 

We next consider the effects of doping on the kagome ice, as shown in Fig. |2](d,e,f) where 
we highlight the vertex configurations in a portion of a system with x = 0.095, with 
the red circles indicating the locations of the doubly occupied traps. For the kagome 
ice, Nl gs and N2 gs denote the two-out, one-in and two-in, one-out ice rule obeying 
states, respectively, NO are the +1 three-out monopoles, and N3 are the -1 three-in 
monopoles. Figure E](d) illustrates the low temperature behavior at T = 1.0. Unlike 
in the square ice, addition of doping sites to the kagome ice has essentially no effect 
on the ground state configuration since it is already highly degenerate. The doubly 
occupied traps effectively inject additional “in”-pointing spins into the system. In an 
undoped sample the number of Nl gs and N2 gs vertices is roughly equal; however, when 
doubly occupied defects are added, the system can absorb the extra “in” spins without 
creating non-ground-state vertices by increasing the fraction of N2 gs vertices in the 
ground state. This flexibility of the ground state is absent in the square ice. As the 
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temperature increases, we observe a continuous increase in the number of N3 monopoles, 
as illustrated in Fig. EX e ) at T = 12; however, there are no compensating NO monopoles. 
Formation of N3 monopoles is energetically favored since there are an excess number of 
N2 gs vertices in the doped ground state. At T = 15, shown in Fig. [2(f), the thermal 
fluctuations are large enough to create both NO and N3 vertices. We find that the doped 
traps in the kagorne ice generally do not act as nucleation sites for additional defects. 
For higher T, the colloids begin to escape from the double-well traps, which sets an 
upper limit on the range of temperature we can study. 

To better understand why the doping has little effect on the kagorne ground state, 
in Fig. [1(a) we plot P G g versus T for kagorne samples with x ranging from x = 0 to 
x = 0.476. In Fig. [1(b), the P G g versus x curves for the kagorne and square ices at T = 0 
indicate that the ground state configuration is lost much more rapidly with increasing 
doping in the square ice than in the kagorne ice. There is a small dip in P G g in Fig. H(a) 
near T = 2.0 for intermediate doping, followed by a local maximum in P G g at T — 5. 
By condensing some of the “in”-pointing spins into N3 vertices, the ratio of Nl gs to N2 gs 
ground state vertices can be shifted back towards its unbiased level. For T < 2.0, the 
decrease in P G g with increasing T in the doped samples occurs as the energy barrier to 
condensation of N3 vertices out of the N2 gs -biased ground state is overcome thermally 
at locations near doping sites where this barrier is suppressed, leading to the formation 
of larger numbers of N3 vertices as the temperature rises. This process halts at T = 2.0, 
the melting temperature of N3 vertices in an undoped system. For 2.0 < T < 2.5, 
the N3 vertices melt back into N2 gs vertices and P G s increases with increasing T. The 
ground state melting temperature is T — 2.5, as shown by the black line in Fig. [4(a) 
for the undoped sample, so for T = 2.5 and above, P G s decreases with increasing T as 
thermally excited defects emerge throughout the sample. 

For T > 14 and low doping, some of the doped states develop higher ground state 
order than the undoped sample. This is more clearly illustrated in Fig. [4(c), where we 
plot P GS versus x at T = 17. Here P G g initially increases with increasing doping before 
reaching a local maximum at x — 0.075 and then decreasing as the doping is further 
increased. This shows that in certain cases the doping can suppress the thermally 
induced creation of NO and N3 defects by locally breaking the degeneracy of the kagorne 
ground state. For comparison, in Fig. [41(d) we plot P G g versus x at T = 17 for the 
square ice, where we find a monotonic decrease in P G g with increasing doping. 

5. Dynamics 

To further quantify the different effects of doping in the square and kagorne ices, we 
examine the hopping rate kh at which the colloids jump from one well to the other 
inside the double-well traps, in units of inverse simulation time steps. We average k^ 
over two colloid populations: k £ is the average hopping rate for all colloids in traps that 
are next to doped, doubly-occupied traps, while k[ is the average hopping rate for the 
colloids in traps that are not next to doping sites. Here, a trap is defined to be next to 
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Figure 5. (a) Hopping rate kh vs T for the square ice system. Blue curves: for 

vertices adjacent to doubly-occupied traps; red curves: k J h for vertices not adjacent 
to doubly-occupied traps. The curves are plotted for x = 0, 0.00945, 0.0236, 0.0472, 
0.0709, 0.0945, and 0.1181, from the bottom red (blue) curve to the top red (blue) 
curve. The hopping rate is enhanced near the doping sites in the square ice. (b) The 
difference in the hopping rates A kh = k l — vs x for T = 8 (black circles), T = 9 
(dark blue up triangles), T = 10 (light green diamonds), T = 11 (dark green squares), 
T = 12 (light blue right triangles), T = 13 (red down triangles), T = 14 (light orange 
left triangles), T = 15 (dark orange up triangles), T = 16 (light purple diamonds), 
and T = 17 (dark purple squares). For small values of x, A kh is negative since the 
hopping rate is highest close to the doping sites, and it approaches zero as the doping 
level or temperature is increased, (c) (blue) and (red) vs T for the kagome ice at 
x = 0, 0.0094, 0.0238, 0.0476, 0.07142, 0.0952, and 0.119 from the bottom red (blue) 
curve to the top red (blue) curve. Here the trend in the hopping rate is reversed, with 
the lowest hopping rates close to the doping sites, (d) Skh vs x for the kagome ice for 
T = 8 to 17 from top to bottom, with the same symbols as in panel (b). 


a doping site if it is part of a vertex that includes the doped trap. In Fig. Efa) we plot 
fcjj and k[ versus T for square ice samples with x ranging from £ = 0 to £ = 0.1181. For 
finite doping, k £ for the traps close to doping sites increases from zero near T — 2, while 
k[ for the traps that are away from doping sites does not begin to increase until T — 5 
to 7. The temperature at which the initial upturn of k[ occurs decreases with increasing 
x, caused when the hopping barrier in an undoped site is depressed by proximity to two 
or more doped sites, an arrangement that becomes more common as x increases. This 
effect appears more clearly in Fig. |6^b) where we plot the spatial distribution of hopping 
rates in a square ice system with x = 0.0709 at T = 12. Here the largest hopping rates 
occur for traps that are immediately adjacent to two or more doping sites. In Fig. [5](b) 
we plot the difference in the hopping rates A kh = k[ — k^ as a function of x for square 
ice samples with T = 8 to 17. For low doping x < 0.2, the magnitude of A kh is largest 
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Figure 6. The spatial distribution of the thermally induced hopping in the entire 
sample at T = 12.0. Black circles indicate the doped traps. Hopping rates are indicated 
on a color scale where red sites have high kh and blue sites have low kh- (a) In the 
kagome ice at x = 0.07142, traps near the doping sites have suppressed hopping rates, 
(b) In the square ice system at x = 0.0709, traps near the doping sites have increased 
hopping rates. 


at low temperatures and gradually diminishes to zero as the temperature increases. For 
x > 0.2 we find a reversal in the relative hopping rates, where the hopping switches from 
being more rapid close to doping sites at low x to being more rapid away from doping 
sites at higher x. This reflects a transition to a regime where the density of doped sites 
becomes high enough to strongly constrain the colloid positions in many undoped sites 
due to random clustering of the doped sites, so that traps close to doped sites have a 
suppressed rate of thermally induced hopping. 

In Fig. [5]^c) we plot the hopping rates k £ and k{ versus T for the kagome ice at 
doping levels x = 0.0095 to 0.119. Here we find the opposite behavior from the square 
ice, with colloids in traps next to doping sites having a significantly lower hopping rate 
than colloids in traps that are away from doping sites. Figure [5jd) shows A kh is largest 
at the lowest temperatures and monotonically decreases both with increasing T and 
with increasing x. In Fig. [6]^a) we illustrate the spatial distribution of the hopping rates 
in the kagome ice for x = 0.0714 at T — 12, where the lowest hopping rates appear 
adjacent to the doped sites. The clustering effect observed for the hopping rate in the 
square ice system, where there is a change in sign of A k^ as the doping level increases, is 
absent in the kagome ice. These results show that doping has opposite effects on square 
and kagome ices. Since the square ice has an ordered ground state, the doping adds 
local frustration which serves as weak spots at which hopping can more readily occur. 
In the kagome ice, the ground state is not ordered, and doping instead lifts the local 
degeneracy to create “hard” spots at which the hopping rate is suppressed. 
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In summary, we have examined the effect of doping on square and kagome artificial spin 
ice systems constructed from colloids in double well traps. Doping is introduced by 
adding an extra colloid to a single trap to create an effective spin that is pointing 
both in and out of the corresponding vertex. For the square ice, doping creates 
additional monopole excitations that serve to screen the doped site from the bulk. As the 
temperature increases, these screening monopoles begin to move away from the doping 
sites by generating additional monopoles or a string of biased ground state vertices. 
The doping sites create local weak spots at which local thermal disordering can readily 
occur as the temperature is increased. For kagome ice, we find the opposite effect, with 
no additional monopoles created by doping since a shift in the population of ground 
state vertices serves to absorb the extra charge introduced by doping, while there is a 
reduction of the hopping rate of colloids in traps adjacent to the doped sites. The doping 
can cause a small increase in the fraction of vertices in the ground state configuration at 
finite temperatures in the kagome ice due to the suppression of the hopping, while for 
the square ice the doping always decreases the fraction of vertices in the ground state. 
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